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ABSTRACT 

We investigate if the discrepancy between estimates of the total baryon mass fraction obtained from 
observations of the cosmic microwave background (CMB) and of galaxy groups/clusters persists when 
a large sample of groups is considered. To this purpose, 91 candidate X-ray groups/poor clusters 
at redshift 0.1 < z < 1 are selected from the COSMOS 2 deg 2 survey, based only on their X-ray 
luminosity and extent. This sample is complemented by 27 nearby clusters with a robust, analogous 
determination of the total and stellar mass inside Rsoo- The total sample of 118 groups and clusters 
with z < 1 spans a range in M500 of ~ 10 13 -10 15 M©. We find that the stellar mass fraction associated 
with galaxies at R500 decreases with increasing total mass as M 50 ° 37±0 04 , independent of redshift. 
Estimating the total gas mass fraction from a recently derived, high quality scaling relation, the total 
baryon mass fraction (/5oo rs+gas = /foo rs + /foo) ^ s f° un d to increase by ~ 25% when M 50 o increases 



from (M) = 5 x 10 3 M© to (M) = 7x 10 14 M©. After consideration of a plausible contribution due 
to intra-cluster light (11-22% of the total stellar mass), and gas depletion through the hierarchical 
assembly process (10% of the gas mass), the estimated values of the total baryon mass fraction are 
still lower than the latest CMB measure of the same quantity (WMAP5), at a significance level of 
3.3ct for groups of (M) = 5x 10 13 M . The discrepancy decreases towards higher total masses, such 
that it is la at (M) = 7 x 10 14 M©. We discuss this result in terms of non-gravitational processes 
such as feedback and filamentary heating. 

Subject headings: galaxies: clusters: general — galaxies: stellar content — cosmological parameters 
— cosmology: observations — X-rays: galaxies: clusters - diffuse radiation 
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1. INTRODUCTION 

The baryon mass fraction is a parameter which can be 
constrained by the primordial light element abundance 
set by the nucleosynthesis at very early times. It has 
been measured to a very high precision from the 5 years 
Wilkinson Microwave Anisotropy Probe (WMAP5) ob- 
servations of the Cosmic Microwave Back ground (CMB), 
givin g a value of /Wmaps = Q,i7i±o,Q09 (|D UI1 klev et all 
2009) 21 . An independent measure of this quantity can 
also be achieved with galaxy clusters. These structures 
are large enough to be representative of the baryon con- 
tent of the universe, which exists mainly in the form of 
X-ray emitting gas and stars. In the absence of dissipa- 
tion, they are expected to provide a baryon mass frac- 
tion ff, comparable to the one m easured from the CMB 
(|White et al.lfl^9l : lEvrard[l99l . 

Galaxy systems appear in a wide range of masses, 
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from ~ 10 to ~ 10 M q. In a hierarchical sce- 
nario dWhite fc FrenidflQQlh the less massive ones, (M< 
10 14 M Q , referred as groups) are the building blocks for 
the most massive ones (clusters). However, the vast ma- 
jority of the attempts to estimate the baryon mass frac- 
tion in ne arby clusters have reported smaller values than 
expec ted (lEttorill2003l:|Lin et al.| [200l iBiviano fc Saluccil 
120061 iMcCarthv et al I I2007T ). In addition this discrep- 
ancy appears to be larger for groups than for clusters 
(|Lin et al.l 2003). Explanations for this invoke physical 
processes which lowe r fb in clusters relat i ve to the univer - 
sal fraction (see e.g. iBialek et all 120011 : iHe et al.l [2006h . 
baryon components that fail de tection by standard X-ray 
and / or optical techniques (see lEttor i 2003; L in &: Mohrl 
l2004f). or a systematic underestimate of f2 m by WMAP 
(IMcCarthv etall[200l 

IMcCarthv et alJ l 20071 ) extensively discuss possible ex- 
planations for the missing baryons. They concluded that 
the observed stellar mass function limits the contribu- 
tion by low mass stars and brown dwarfs to a negligible 
fraction of the total stellar mass; furthermore they rule 
out a contribution by large amounts of centrally concen- 
trated gas, on the bases of inconsistencies with current 
X-ray data and the assumption of hydrostatic equilib- 
rium. Consideration of the so called intra-cluster light 
(ICL) results into a discrepancy at the 3.2cr level with re- 
spect to WMAP3 acro ss the mass range 6x 10 13 -10 15 M Q 
((Gonzalez et al.|[2007t ). As discussed by these authors, 
systematics may help reconciling their results with the 
WMAP estimate. 

In this respect, the correct determination of the 
gas mass fraction may be crucial. In fact, studies 
of the individual baryon components (stars associated 
with galaxies and gas) have shown that the stellar 
(/focT = ^50o rs /^50o) and gas mass fractions within 
R-500 22 (/foo = Mfoo/M5oo) exhibit opposite behaviours 
as a function of the total system mass. In particular 
clusters have a highe r gas mass fraction t h an groups 
(IVikhlinin et alj 120061; I Arnaud et all l2007t iSun et al " 
2009), but a lower stellar mass fraction dLin et al 



2003). This has been interpreted as a difference in the 



star formation efficiency between groups and cluste rs 
(|David et al.lll990t iLin et al.ll200l iLagana et al"1l2008f ) 

On the other hand the mass dependence of the gas 
fraction and the discrepancy between the baryon mass 
fraction in groups/clusters and the WMAP value can be 
understood in terms of non-gravitational processes. In 
fact AGN-heating (which can drive the gas outside the 
potential well) or gas pre-heating (which inhibites the 
gas from falling towards the center of the potential) can 
explain the lack of gas within rsoo in groups. Therefore 
groups appear as the critical systems to assess the univer- 
sality of the baryon fraction, and to understand complex 
physical processes affecting both the gas and the stellar 
components. 

Little work has been conducted on estimation of the 
baryon mass fraction at the group regime, mainly be- 
cause of the lack of groups in existing catalogues and the 
difficulty of estimating masses for the individual compo- 

22 R A (A=500, 200, 2500) is the radius within which the mass 
density of a group/cluster is equal to A times the critical density 
(p c ) of the Universe. Correspondingly, Ma = A p c (z) (4 7r/3)i? A 
is the mass inside Ra- 



nents and the total. An insufficient sampling of the range 
in total mass spanned by groups and clusters is problem- 
atic for studying their overall properties in terms of mean 
and scatter of the population 23 . A galaxy group/cluster 
is the result of the assembly history of the dark matter 
halo, as well as of the star formation processes affecting 
the gas. Both processes lead to multivariate outcomes 
and produce a large intrinsic scatter in the distribution 
of the observed properties of groups and clusters. There- 
fore it is essential to have a large enough sample to be 
representative of the population, and unbiased by selec- 
tion effects, to be able to investigate the mean trend 
precisely. 

Once such a sample is available, interesting questions 
to address are: (1) How does the stellar mass fraction 
behave across the total range of masses? (2) Does the 
relation between the stellar mass fraction and the to- 
tal system mass evolve with redshift? (3) How does 
the gas mass fraction change as a function of the sys- 
tem total mass? (4) Is the total baryonic fraction in 
groups/clusters of galaxies consistent with the WMAP5 
value? 

In this paper we select the currently largest X-ray se- 
lected sample of groups from the COSMOS 2 deg 2 survey 
which consists of 91 high-quality systems at 0.1 < z < 1. 
Existing observations currently do not give constraints 
on the evolution of the baryonic components in individ- 
ual systems at z>0.1. Our data allow us to put con- 
straints on the redshift evolution of the average stellar 
fraction with mass, which we find to be consistent with 
zero (§4.2). Observational constraints on the evolution 
of the average gas mass fr action also sugge st zero evolu- 
tion in the cluster regime ([Alien et al.l |2004t). We assume 
that this is applicable to our groups in the absence of ob- 
servations to the contrar y and we note that s imulations 
support this hypothesis (|Kravtsov et alJ l2005). 

We comple ment our sample with 27 nearby clusters in- 
vestigated bv lLin et al.l (|2003l ) in order to achieve a span 
of two orders of magnitude in total mass (10 13 <M< 
10 15 M ). In <j3]the total mass of stars associated with 
galaxies is directly determined for each group, and we 
investigate the relation between the stellar mass fraction 
and the total mass of the system. In §4 we combine the 
stellar mass fraction estimates with the most recent de- 
termination of the relation between gas mass fraction and 
total mass based on a compilation of 41 local (z <0.2) 
X-ray groups and clusters , span ning the same range in 
mass as ours (jPratt et al.l [2009) , and we compute the 
total baryon fraction. We discuss results in §l5j 

We adopt a ACDM cosmological model (f2 m = 0.258, 

O a = 0.742) wit h H n = 72 km s" 1 Mp c " 1 , consistently 
with WMAP 5 (jDunklev et alj l2009t iKomatsu et all 
2009). Unless otherwise stated all quantities are esti- 
mated at an overdensity of 500. 

2. THE SAMPLE 
2.1. The COSMOS survey of groups /poor clusters 
The Cosmic Evolution Survey (COSMOS, 



iScoville et all l2007ah was designed to probe how 
galaxies, active galactic nuclei (AGN), and dark mat- 



23 The conclusions of ILin et all l|2003T ) and IGonzalez et al.l 
1)20071 ) . for example, are based only on, respectively, 27 and 23 sys- 
tems, but only 3 and 5 of them are less massive than 10 14 h~ 1 Mq. 
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ter evolve together within the large-scale structure. 
The survey is based on multi-wavelength imaging 
and spectroscopy from X-ray to radio wavelengths 
and covers a 2 deg 2 area, includ i ng H ST imaging of 
the entire field (|Koekemoer et all [2007). Large-scale 
structures in the COSMOS field have been characterized 
in ter ms of galaxy overden sity using photometric red- 
shifts (IScoville et al.1 l2007bfl . weak lensing convergence 
maps dMassev et all 120071 ) . diffuse X-ray emission 
Finoguenov et al. |2007[ ) and a combination of these 



Guzzo et al.ll2007D . In particular, the entire COSMOS 
region was imaged t hrough 54 overlapping XMM- Newton 
pointings (1.5 Ms, iHasinger et al. 2007). Additional 
Chandra observations (1.8 Ms. lElvis et al. l2006f ) mapped 



the central region to higher resolution. 

In this study we use X-ray detection, gravitational 
lensing signal, optical photometric and spectroscopic 
data of the clusters and groups identified in the COS- 
MOS survey . The X-ray dat a redu ction is described 
in detail in iFinoguenov et all (|2007l ) and Finoguenov 
et al. (in preparation). From a composite mosaic of 
the XMM-Newton and Chandra X-ray data, it has been 
possible to detect and measure the flux of extended 
sources (i.e., groups and clusters) down to a limit of 
10~ 15 erg s _1 cm~ 2 , as described in the corresponding 
catalogue (Finoguenov et al. in preparation). Extended 
source detection was based on a wavelet scale-wise re- 
construction of the image, as described in Vikhlinin et 
al. (1998b), employing angular scales from 8" to 2.1'. 
Clusters and groups of galaxies were effectively selected 
by the spatial extent of their X-ray emission, following 
the approach of Rosati et al. (1998), Vikhlinin et al. 
(1998b), and Moretti et al. (2004). The cluster detection 
algorithm consists of: (1) removal of the background, (2) 
detection of AGN, (3) removal of AGN flux from large 
scale, and (4) search for extended emission. As a result, 
a total of 219 X-ray extended sources were identified in 
the redshift range < z < 1.6; they span the rest-frame 
0.1-2.4 keV luminosity range 10 41 < L x < 10 44 erg s -1 , 
which is typically populated by groups and poor clusters. 

Quality flags tag individual systems. Flag 1 is assigned 
to objects whose center corresponds to the X- ray peak 
of the source, while flag 2 objects have their spectral 
extraction region redefined to include only their robust 
association with a unique optical system. A redshift was 
assigned to each candidate X-ray group/cluster, corre- 
sponding to the mean of the photometric redshift (photo- 
z) distribution of the red-sequence galaxies as identified 
in Tanaka et al. (in preparation), if present, and ly- 
ing within the X-ray overdensity contour region. This 
redshift is checked against the available spectroscopic 
redshifts most ly provided by the zCOSMOS spectro- 
scopic survey (iLillv et alJl200l . The presence of a red- 
sequence is not required for the group/cluster detection: 
if no overdensity of red sequence galaxies is found in the 
photo-z space, the spectroscopic data only are checked 
for the presence of a galaxy overdensity in the same area. 
Flag 3 is assigned to high-z (z>l) not spectroscopically 
confirmed candidate groups. Flag 4 is assigned when 
multiple optical counterparts are present within the X- 
ray overdensity contour region. In this study only sys- 
tems with quality flag 1 or 2 are considered. 

The galaxy-group detection is irrespective of any op- 
tical characteristic, being based only on the presence of 



an X-ray extended source. The X-ray selection is an ap- 
proximate selection by hal o mass, due to the tight X-ray 
luminosity-mass relation ([Pratt et al.ll2009T ): in this re- 
gard our selection is thus unbiased with respect to both 
the optical properties of the groups in our sample and 
the X-ray characteristic of the systems. 
The purposes of the present study lead us to intro- 
duce three further selection criteria: (1) only candidate 
groups/clusters detected in X-rays with a significance 
higher than 3cr on the flux determination are considered. 
Selection of the most robust candidates minimizes con- 
tamination by loose galaxy aggregations or superposi- 
tion of AGN along the line of sight. (2) Only X-ray 
extended sources with Lx > 10 42 erg s _1 are considered, 
in order to limit co ntamination from starburst galaxies 
([Grimm et al.l I2003D or field ell iptical galaxies with X- 
rav halos (|Diehl fc Statierll2007h . (3) We limit the red- 
shift range to 0.1 < z < 1.0 , where photo- z o f individ- 
ual galaxies are most robust ([Ilbert et al.ll2009l ): further- 
more, in this range the quality of the photo-z is equiva- 
lent to that of low resolution spectroscopy. 

Figure [T] reproduces the X-ray luminosity distribu- 
tion as a function of redshift for the candidate X-ray 
groups/clusters within z= 1 (151 out of 219 systems). 
The flag 1+2 sample selected for this s tudy contains 114 
object s, of which 44 were present in IFinoguenov et aTl 
(2007). It contains only 3 systems at z < 0.2 (Fig- 
ure [l}, which is the redshift range co vered by analo- 
gous studies on ,fb in groups/clusters ([Lin et al.1 120031 : 
iGonzalez et al.l2007t l. On the other hand, it contains sys- 
tems with particularly low X-ray luminosities (i.e., with 
10 42 < L x < 5 x 10 42 erg s" 1 ), though only for z < 0.5. 
The sample considered in this study is reduced to 91 ob- 
jects after removal of 23 groups with unreliable estimates 
of the total stellar mass in galaxies f§ 13. 2|) . Out of these 
91 candidate groups/poor clusters, 51 are already spec- 
troscopically confirmed (i.e. are associated with at least 
3 galaxies with similar spectroscopic redshifts). 

2.2. COSMOS X-ray-selected groups/poor clusters: 
total mass estimate 

In the original catalogue (Finoguenov et al. in prepa- 
ration), M200 is computed using an LX-M200 relation es- 
tablished via the weak lensing analysis in Leauthaud et 
al. (in preparation). Briefly, the COSMOS group sam- 
ple is divided into nine bins that span the redshift range 
0.1 <z< 0.9 and with 10 418 <L X /E(z)< 10 43 5 erg s" 1 , 

where the function E(z)= y/Sl m (l + z) 3 + H,\ represents 
the Hubble parameter evolution for a flat metric. Only 
systems with a clear, visually identified BCG are used for 
this analysis, to minimize issues due to incorrect center- 
ing. For each bin, the weak lensing signal is calculated 
from r ~ 50 kpc to r ~ 3 Mpc in logarithmically spaced 
radial bins. A weak lensing signal is detected all the 
way to 3 Mpc ensuring that the lens density is probed 
well beyond the virial radius. The results are fit with 
a p arametric model whi ch is the sum of a NFW pro- 
file ([Navarro et al.l H997rj and a point-source term due 
to the mass of the central BCG. The theoretical relation 
between mass and concentration from IZhao et all ([2008) 
has been used in the fit for the NFW component and the 
mean stellar mass of the central BCG's is used in order to 
scale the point source term. A comparison between the 
relation obtained from the combination of the the COS- 
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Fig. 1. — Rest-frame 0.1—2.4 keV luminosity vs. redshift for the 
151 COSMOS candidate X-ray groups/clusters at 0.1< z <1.0. 
Filled circles mark the 91 objects considered in this study: dark or 
light grey identifies objects with flag 1 (45) or 2 (46), respectively. 

MOS data and cluster data from Hoekstra et al . 2007 
is consistent with that obtained by Rykoff et al. (2008) 
based on SDSS data. We adopt the following functional 
form for the Lx~M relation, 



Mwo E(z) 
M 



= A 



L x E(z)- 



L 



XX) 



(1) 



where M = 10 13 7 M , L x ,o = 10 42 7 erg s" 1 . Fitting 
only the COSMOS data yields the best fit parameters 
log 10 (A) = 0.106±0.053 and a = 0.56±0.12 (cited errors 
are statistical only). Further details regarding the weak 
lensing analysis in COSMOS can be found in Leauthaud 
et al. in preparation. 

The baryon fraction in groups/clusters can be stud- 
ied at any radius, though it is desirable to study it at 
the largest radius possible with respect to the virial ra- 
dius of the system because of the radial dependencies of 
the different components. The largest radius for which 
relia ble X-ray hydrosta t ic masses are available is Rfsno 
(e.g.lArnaud et al.ll2005l : IVikhlinin et~aLll2006t ISun et all 
2009) . Hereafter we use M500 instead of M200 > to enable 
a comparison at the same radius with previous studies 
on /b in nearby groups/clusters. The catalogue value of 
M 2 oo is converted into M 50 o assuming an NFW profile 
with a constant concentration parameter (c =5). 

2.3. COSMOS galaxies: multiwavelength photometry 
and photometric redshifts 

The COSMOS area ha s been imaged in 30 ba nds in- 
cluding broad- fSUBARU lTaniguchi et all2007at CFHT 
McCracken et al. in preparation), medium-, and narrow- 
bands (SUBARU; Taniguchi et al. i n preparation), 
rangi ng from the far-ultraviolet (G ALEX Zamoisk i et al.l 
l2007f ) to the mid-infrared (gpifeer ISanders et al.ll2007f ). 
This multiwavelength dataset is collected in a master 
photometric catalogue. Capak et al. (in preparation) 



discuss in detail source detection and extraction of pho- 
tometry. The COSMOS photometric catalogue is com- 
plete down t o a to tal i-b and magnitud e of 26 .5 AB mag. 
Illbert et al.l (|2009f) and ISalvato et ail (|2009D computed 
highly reliable photometric redshifts with unprecedented 
accuracy for a survey this large, owing to the extraor- 
dinarily large number of photometric bands. Redshifts 
were attributed to individual galaxies v ia a standard x 2 
fitting procedure (jArnouts et alJ 120021 ) encoded in Le 
Phare 24 , written by S. Arnouts and O. Ilbert. Best- 
fit solutions from this photo-z algorithm were trained 
on a composite spectroscopic sample of objects brighter 
than zab = 25 (see table 3 in Illbert et al.ll2009l ). mostly 
made of ~4,000 bright galaxies (i.e., with ^ab < 22.5) 
observed as part of the zCOSMOS spectroscopic survey 
(jLillv et al.ll2007h . Comparison of photometric and spec- 
troscopic redshifts gives a typical r.m.s. scatter of the 
photo-z's equa l to crphoto-z = 0-0 2 x (1 + z) for zab < 25 
and z < 1.25 (jllbert et al. 2009). In the presence of X- 
ray emitting objects (AGNs ), photometric redshi fts were 
independently estimated by ISalvato et alJ (|2009f h 

As a by-product of the photo-z determination, spec- 
troscopic types were attributed to individual galaxies on 
the basis of their best-fit broad-band spectral energy dis- 
tributions (SEDs). This information is used to estimate 
the stellar mass of a galaxy, which is obtained from the 
conversion of the Ks-band luminosity (jllbert et al.ll2009h 
using an evolving galaxy-type depe ndent stellar mass-to - 
Ks-band luminosity ratio M/Lk s (jArnouts et al.l 1 2007). 
This relation h as been establi shed using a Salpeter initial 
mass function (|Sah>eterlll955l ) . Stellar masses of individ- 
ual galaxies are contained in the COSMOS photometric 
catalogue; the fractional error on the stellar mass of a 
galaxy is typically equal to 3 4% , and is domina ted by 
the mean scatter on M/Lks (jArnouts et al.|[2~007f ). 

This uncertainty pertains to the aforementioned 
method of estimating stellar masses. Individual galaxy 
stellar masses may differ by a factor 2-3, depend- 
ing on the method used to estimate the mass (e. g 
iLonghetti fc Saraccol [20091 : iKiipcu Yoldas et al.l l200l . 
This uncertainty is the product of several factors; it 
mostly reflects the range of assumptions in differing mod- 
els as for the star-formation history (e.g., single burst 
vs. multiple bursts vs. continuum star-formation ac- 
tivity) and the attenuation of stellar light by dust (e.g., 
starburst-like vs. normal star-forming disc-like). In ad- 
dition, it results from different implementations of com- 
plex physics, such as the asymptotic-giant-branch phase 
of stellar evolution and metal enrichment). This scatter 
does not reflect the uncertainty of the present method, 
which is 34% for individual galaxies as detailed above. 
This latter value is the uncertainty we attribute to indi- 
vidual galaxy stellar masses in the present study. 

2.4. Nearby clusters 

The COSMOS sample is mostly composed of groups. 
Therefore we complement it with a sample of 27 nearby 
X-ray select ed clusters wit h sufficiently deep 2MASS 
photometry (|Lin et al.|[200l LMS03) to estimate accu- 
rate stellar masses. The total and stellar masses were 
derived by LMS03 in a manner consistent with ours. In 
particular, the total cluster mass is estimated from an 

24 www.lam.oamp.fr/arnouts/LE_PHARE.html 
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M500-TX relation. The stellar masses are estimated from 
the total K band luminosity of each cluster, assuming an 
average stellar mass-to-light ratio which takes into ac- 
count the varying spiral galaxy fraction as a function of 
the cluster temperature. 

LMS03 provide estimates of the total gas fraction ob- 
tained from either X-ray data or from a scaling rela- 
tion; we use inste ad the most recent scaling relations of 
iPratt et~al . (2009), based on hydrostatic mass estimates, 
in order to reduce systematic effects. We apply this both 
to our sample and the one of LMS03. 

3. DATA ANALYSIS 

3.1. Galaxy stellar mass function: completeness and 
extrapolation 

The low-mass end of the galaxy stellar mass function of 
the individual COSMOS groups/poor clusters is probed 
to different extents by observations, since these systems 
span a rather large redshift range (0.1 < z < 1). In order 
to achieve a common footing, the completeness in galaxy 
absolute magnitude (stellar mass) of the sample must be 
understood. 

First, we divide the sample into two redshift bins (0.1- 
0.5 and 0.5-1.0) containing a similar number of objects, 
since the cosmic stellar mass density is ob served to drop 
by a factor of 2 from z ~0 to 1 in the field (|Wilkins et al.l 
2008 and references therein). The completeness mass is 
estimated at z=0.5 and z=1.0 from a fit of its behaviour 
as a function of redshift, obtained using a sampling of 0.1 
in redshift as follows (Bolzonclla et al. in preparation). 
Firstly we derive the stellar mass (M\i m ) that each object 
would have if its apparent magnitude was equal to the 
sample limit magnitude (i.e. iAB=25), viz. , 



log Af lim = log A/ + 0.4 x (i AB - 25.0), 



(2) 



where M is the stellar mass of a galaxy with apparent 
magnitude zab- Secondly we derive the 95% percentile 
of the distribution in Mi; m for galaxies in the lower 20% 
percentile in magnitude (i.e. iAB > 23.6) in each bin of 
0.1 in redshift. Finally a fit to the corresponding envelope 
as a function of redshift is performed for 0.1< z <1.0; the 
ensuing values represent the stellar mass completeness as 
a function of redshift for our sample. Figure [2] illustrates 
the behaviour of the stellar mass completeness as a func- 
tion of redshift. 

For instance, the stellar mass completeness at z=l 
(M comp i = 1O 1O - 4 M ) is about an order of magnitude 
lowe r than the so-called "transition" stellar m ass at z < 1 
(e.g. lBundv et al.ll2005t fPannella et al.ll2006D . This con- 
firms that a rich mixt ure of morpho logies and, thus, 
star-formation histories (Sandagc 1986) is present among 
the member galaxies of the COSMOS X-ray selected 
groups/poor clusters. 

We compute the total stellar mass associated with 
galaxies of a given system as follows. We first add the 
stellar masses of galaxies more massive than the com- 
pleteness mass (at z=0.5 or 1) for which membership to 
a given group/poor cluster is determined (as described 
in H3.2.1[) . Taking into account the mass of the individ- 
ua l galaxies, rathe r than their statistical distribution (as 
in iLin et al.ll2003T ). becomes increasingly important for 
groups, where the BCG is a large fraction of the total 
stellar mass. 




0.6 
redshift 

Fig. 2. — The completeness stellar mass for our sample is com- 
puted from the fit (black dashed line) to the 95% percentile of 
the distribution in Mu m (see text) for galaxies in the 20% lower 
percentile in magnitude (grey circles) as a function of redshift. 
The black dots represent the stellar masses for all galaxies with 
i\B < 25. To reduce the plot size, we plot only one point in ten. 

The contribution from less massive galaxies is esti- 
mated in a statistical manner from the composite stellar 
mass function (Giodini et al. in preparation), which can 
be robustly obtained only within two broad redshift bins 
(0.1< z <0.5 and 0.5< z <1.0). The stacked stellar 
mass function for systems falling in ea ch redshift bin i s 
fitted with a single Schechter function (|Schechterlll976h : 
the correction factor for stellar masses lower than the 
completeness mass, down to ~10 8 M Q (typical mass of a 
dwarf galaxy), is given by: 



/m° J{M)-MdM 



C 3 f(M) 



AM 



(3) 



where M comp i is the completeness mass for the given red- 
shift range. The fractional contribution to the total stel- 
lar mass budget of galaxies with 10 8 M < M< M comp i 
corresponds to -9% (-1%) at redshifts 0.5-1.0 (0 . 1-0.5) . 
These values are almost negligible, as in lLin et~a l. ( 2003), 
which confirms that the total stellar mass associated with 
galaxies can be achieved almost directly from the data 
for our sample of X-ray-selected groups/poor clusters at 
0.1 < z < 1. 

3.2. Total stellar mass (in galaxies) 
3.2.1. Statistical membership 

As a first step, we estimate a projected total stellar 
mass, which is the sum of the stellar masses of all poten- 
tial member galaxies down to the completeness mass of 
either redshift bin to which a group belongs (i.e. 0.1-0.5 
or 0.5-1.0). Candidate members are defined as all the 
galaxies within a projected distance equal to R500 from 
the X-ray centroid of a group/poor cluster and within 
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0.02x(l+z) from its redshift (given in the X-ray cata- 
logue). Then we perform a foreground/background cor- 
rection by measuring the total stellar mass of galaxies 
contained in 20 circular areas which have the same ra- 
dius as R500 and have photometric redshifts consistent 
with that of a given system within the errors. These ar- 
eas do not overlap either with the group or with other 
groups at the same redshift and are chosen to represent 
the coeval field environment. Field galaxies are selected 
in redshift and stellar mass following the same criteria as 
for the selection of potential member galaxies previously 
described. The mean and the standard deviation of the 
distribution of the total stellar masses computed in the 
20 regions are taken as the value of the stellar mass as- 
sociated with the foreground/background and its uncer- 
tainty, respectively. Finally, the foreground/background 
value is subtracted from the initial estimate of the total 
stellar mass of the system. 

If the error on the foreground/background value is larger 
than half of the estimated total stellar mass content in 
galaxies of a given system, this system is removed from 
the sample. Obviously a system is excluded also if the 
foreground/background correction exceeds the estimated 
total stellar mass content in galaxies. The variance on 
the total stellar mass budget in galaxies for a system is 
given by the sum in quadrature of the background un- 
certainty and the error on the total stellar mass of the 
galaxies of the system. 

Furthermore we checked the influence of masked areas 
on the reliability of the computed total stellar masses 
of individual groups. A region of the COSMOS area is 
masked when the image quality is poor owing to differ- 
ent reasons (e.g. field boundary, saturated stars, satellite 
tracks and image defects). For galaxies with elliptical- 
like SEDs reliable photo-z's can be determined also in 
masked areas; therefore early-type galaxies falling in 
masked areas are considered. On average, the contri- 
bution of these objects to the stellar mass budget of a 
group is not expected to be negligible. In fact, in 30 out 
of 37 cases where early-type galaxies falling in masked 
areas are retrieved, the new stellar mass fraction (com- 
puted in £)4.1I ) is consistent with that of other groups 
with the same M500, whatever the redshift. Conversely, 
late-type galaxies falling in masked areas are not consid- 
ered and the impact of this choice is tested a posteriori. 
For 23 out of 114 groups the number of statistically es- 
tablished member galaxies is less than 6 and the total 
stellar mass is systematically lower than the mean for 
groups of similar total masses, irrespective of M 50 o 25 ■ 
These 23 objects span the entire total mass range and 
their exclusion does not affect our results on the stellar 
mass fraction; at the same time, the scatter in the stellar 
mass fraction decreases by 30% 26 . Only the resulting 
sample of 91 galaxy systems with at least 6 members, 

25 This tells that 5 members only is insufficient to determine 
the stellar mass budget of a group. In fact, when the total stellar 
mass or luminosity of a system is computed from a population of 
discrete sources, the scatter in the ensuing value turns out to be 
non linear when the number of discrete so urces becomes small ( e.g. 
of order ten or less), as demonstrated by Gilfanov ct al. (2004) in 
an analogous application. 

26 Nevertheless these objects are potentially an interesting sub- 
population characterized by an extremely slow build-up of stellar 
mass. Further optical follow up will help to better assess their 
properties. 



spanning two orders of magnitude in X-ray luminosity, 
is considered in the following analysis; hereafter it is des- 
ignated the COSMOS X-ray selected group sample. 

3.2.2. Deprojection 

The total stellar mass in galaxies so far estimated refers 
to a cylindrical section of the system projected onto the 
plane perpendicular to the line of sight. We therefore 
need to deproject the total stellar mass from two to 
three dimensions. The average galaxy distribution is de- 
scribed by a projected NFW profile in two dimensions 
(jBartelma nn 1996: iNavarro et al.lfl997h : 



E(aO 



2/O s r s 



/(*) 



(4) 



where 



1 arctan J ^rfj (x > 1) 



1 Tx^i 

1 - arctanh J (x < 1) 

{x = l) 



(5) 



and as a generalized NFW profile in three dimensions 



p(x) 



x(l + x) 2 



(6) 



In both equations the radial coordinate x is the radius 
in units of a scale radius r s , x = r/r s . The scale radius 
corresponds to the ratio between R200 an d the concen- 
tration parameter c for the system. An average profile is 
produced using all 91 systems in our final sample, with a 
central density normalized to the number of groups. This 
high signal-to-noise, average two-dimensional galaxy dis- 
tribution is best-fitted by a two-dimensional NFW profile 
where r s = 0.27 i?200- The average radial profile is shown 
in Figure [3] together with its best fit (with a reduced 
X 2 value equal to 1.2). We remark that our aim is not 
to compute the concentration parameter of the galaxy 
distribution for individual systems, otherwise we should 
take into account the scatter in the evolution of the con- 
centration parameter as a function of redshift. Instead 
we want to compute an average correction for projection 
of the mass profile of a system as calculated in £13.2.11 

Using the best-fit values, we compute correction factors 
by integrating the average profile out to R500 : 



dpf = 



J Q H 500 p(r) ■ 4irr 2 dr 



f* 500 E{r)-2xrdr 



(7) 



The deprojected total stellar mass of a system is then 
given by 

Mf*™ = dpf x M$£ m , (8) 

where dpf — 0.86 is the correction factor. 

4. RESULTS 

4.1. Stellar mass budget (galaxy component) 

Figure 2] shows the behaviour of the total (deprojected) 
stellar mass in galaxies within R500, M|QQ rs , as a func- 
tion of the total mass M 500 for the 91 COSMOS X-ray 
selected groups. The distribution in Figure H] exhibits 
a rather well defined trend, although a large scatter is 
present, especially at low masses, where values can range 
by a factor of 10 at a fixed total mass. Part of this large 



Stellar and total baryon mass fraction in groups and clusters since z=l 



7 



100 



A 

w 

V 



10 



K 



ft, 



0.1 



1.0 



R/R 



200 



Fig. 3. — Radial profile of the average number galaxy density 
for the 91 COSMOS groups/poor clusters. The dashed line shows 
the best fit NFW profile (c~4). The unit of the surface density 
is number per area in unit of ttRjqq and normalized to the total 
number of systems. 

scatter may have a physical origin: different merging his- 
tories produce different to tal mass-to-light r atios for fixed 
total assembled mass (cf. lSales et all 12001 . 

We fit the relation between total stellar mass in galax- 
ies and total mass for all 91 systems and for the 45 flag=l 
groups only. Since the distribution in Figure 0] exhibits 
an intrinsic scatter larger than the errors on the individ- 
ual points, the fit is performed using the weighted least 
sq uare with intrinsic scatter (WLSS) method discussed 
in lPratt et al.l (|2006f ). This algorithm takes into account 
uncertainties on both stellar mass and total mass and 
the presence of intrinsic scatter in the data. There is a 
robust correlation between M|gQ rs and M 500 in the COS- 
MOS X-ray selected groups: 



ji ,f stars 
JW 500 



(0.30 ±0.02) x 



M: 



113 K2 1 



0) 



,5 x 10 1 

where a=0.81±0.11 for the entire sample and 
a=0.72±0.13 for the flag=l subsample, and the (loga- 
rithmic) intrinsic scatter is equal to 35% in both cases 27 . 

Fitting the stellar-to-total mass ratio vs. total mass of 
the system for the full sample of COSMOS X-ray selected 
groups only we find 



JUT = 5.0t°;J x 10- 



AL 



.-,1111 



5 x 10 13 M, 







-0.26±0.09 



(10) 



A fit to the Flag=l sample gives equivalent results. Re- 
markably the relation between the mass fraction of stars 

27 This result is robust against the presence of a pair of groups 
which are detected at the same redshift, but with a separation of 
the order of Rsoo- The two objects of this pair lie above the best-fit 
relation reproduced in Figure 2] perhaps as an effect of a bias in 
their estimated total stellar masses in galaxies. However, new fits 
performed after excluding these two groups give the same results 
as the previous ones. 
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Fig. 4. — Total stellar mass in galaxies vs. total mass for the 
91 COSMOS X-ray selected groups/poor clusters. Filled (empty) 
grey circles identify objects with fiag=l (2). The dashed (dotted) 
line represents the best fit relation derived for fiag=l (all) groups 
(see equation [9} derived taking into account uncertainties in both 
quantities and the intrinsic scatter of the relation. 



TABLE 1 

The best fit parameters for the relation 
between stellar mass fraction and total 
mass (eq. [lo] and eq. [tt)l for three samples 
considered. data were fitted with a 

POWER LAW /|^ rs = N(M 500 /5 X 10 13 M Q ) a . 



Sample 


Log(N) a 


slope 


COSMOS flag=l 


-1.35±0.01 


-0.33±0.12 


COSMOS flag=l+2 


-1.35±0.01 


-0.26±0.09 


COSMOS+LM03 


-1.37±0.01 


-0.37±0.04 



in galaxies and the total mass of the system for the COS- 
MOS X-ray selected groups is consistent within the er- 
rors with the o ne fou nd in nearby clusters by LMS03 and 
lLagana et al.l (|2008l ). We now extend the range of total 
masses using the results from local clusters selected by 
LMS03, converting their measurements to our cosmol- 
ogy. Since these authors do not give the uncertainties 
associated with their total mass estimates, we assign a 
fixed fractional total mass uncertainty equivalent to the 
mean of that for the COSMOS groups (-30%). The best 
fit of the combined sample is 

„, \ -0.37±0.04 

M500 



j-stars _ c n+0.1 i r>- 
/500 — °- U -0.1 X iU 



5 x 10 13 Mi 



(11) 



with a typical logarithmic intrinsic scatter of —50%. The 
data and best fit relations are shown in FigureO To 
better elucidate trends with total mass, we divided the 
data set into five logarithmic bins of equal size in total 
mass, and computed the mean and standard deviation 
of the values of the mass fraction of stars in galaxies 
in each bin using the biweight estimators of iBeers et al.l 
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Fig. 5. — Left panel: stellar-to-total mass ratio vs. total mass for the combined sample of 91 COSMOS X-ray selected groups (same 
symbols as in Figure [4j plus 27 nearby clusters of LMS03 (empty squares). The dashed line represents the best— fit relation derived for 
fiag=l groups of the COSMOS sample and the dotted line represents the fit to all COSMOS groups. The solid line shows the best fit 
relation for all COSMOS groups plus local clusters. All fits are derived taking into account uncertainties in both quantities and the intrinsic 
scatter in the relation. The ensuing fit parameters are given in Table [T] The large points with error bars show the biweight mean and 
standard deviation of these data binned in 5 loga rithmic bins in total mass. Right panel: gas f raction a s a function of t he system mass 
from a combined sample of 41 clusters and groups (IVikhlinin et al.ll200l V06; lArnaud et al.ll2007l, APP07; [5un et al.|[2009l, S08). The solid 
line is the best fit relation /|qq 
bins of total mass. 



M|?qq . The large points with error bars show the mean and standard deviation of these data binned in 5 



TABLE 2 

Measured values for /|qq and /§Q rs as in Figure [6] 
Uncertainties correspond to the standard deviation of 
the mean (see text for details). 



M 5 oo/ [h^Me] 


pstars 
: 500 


rgas 
r 500 


rStars+gas 
'500 


2.1e+13 


0.062±0.005 


0.074±0.028 


0.136±0.028 


5.1e+13 


0.045±0.002 


0.068± 0.005 


0.113±0.005 


1.2e+14 


0.036±0.004 


0.080 ± 0.003 


0.116±0.005 


3.0e+14 


0.021±0.002 


0.103 ±0.008 


0.124±0.009 


7.1e+14 


0.019±0.002 


0.123 ±0.007 


0.141±0.007 



(1990); they are relatively large, which gives a measure 
of the heterogeneity of the population. The large points 
with error bars show the trend of these binned data with 
total mass: there is good agreement with the best fitting 
regression line to the unbinned points, as expected. 

4.2. Evolutionary considerations 

Finally we inspect the presence of evolution of the rela- 
tion between /foo rs and M 50 o by considering only systems 
at z<0.5 (we cannot fit the relation for the high redshift 
systems since they do not cover a sufficient range in to- 
tal mass). The ensuing fit is fully consistent with that 
obtained for the entire sample within the uncertainties. 

We can put a constraint on the possible evolution of 
the relation by evaluating the change in the mean of fstar 
for massive systems (M500 > 5 x 10 14 Mq) in two redshift 
bins (z<0.5 and z>0.5). The average f s t ar changes from 
0.031±0.013 at z<0.5 to 0.039±0.019 at z>0.5, a less 
than one sigma difference in mean values. Even taking 
the maximum distances between the two values given the 



uncertainties, the stellar mass fraction does not change 
by more than 35%. 

A second way to investigate a possible evolution of the 
stellar mass fraction in galaxies is to plot the ratio of 
the stellar fraction to the mean relation as a function of 
redshift (r/(z) = f s tar(z)/ (fstar))- Using the same five 
bins in total mass as above, no trend in r/(z) is evident. 
However a fit of r/(z) gives a robust upper limit on the 
evolution over the maximum redshift range (0-1) of 40%. 
Taking the median redshift of each redshift bin (0.22, 
0.72), the upper limit on the evolution of the stellar frac- 
tion is less than 20%. This number is consistent with the 
upper limit on the evolution of the relation between to - 
tal star fraction and M500 given bv lBalogh etai] ((2008). 
Therefore we conclude that our data do not support the 
existence of a significant evolution in the zero-point and 
slope of the /|Qo rs -M 50 o relation between redshifts and 
1. 



4.3. The total baryon mass fraction 
4.3.1. The gas mass fraction 

In order to determine the total baryon mass fraction in 
individual systems, we need to estimate the amount of 
baryons in the form of hot gas which make the intra- 
cluster medium (ICM). Unfortunately, this cannot be 
achieved from most of the existing X-ray observations of 
the total sample because their signal-to-noise is insuf- 
ficient for the purpose. Therefore, we have to resort to 
an estimate of the mean trend of the gas mass fraction 
as a function of M 50 o established from an independent 
sample of well observed gro ups and clusters at z < 0.2, 
selected from the samples of lVikhlinin et aLll2006l (V06), 
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lArnaud et all [20071 (APP07) and lSun et all [20091 (S08). 
These authors computed gas mass fractions at R500 from 
hydrostatic mass estimates for 10 (V06), 10 (APP07) and 
21 (S08, including the best quality tiers 1 and 2 systems) 
clusters and groups, respectively. The combined sam- 
ple contains 41 systems and spans the total mass range 
1.5xl0 13 -l.lxl0 Mq. After conversion to a common 
cosmology, a fit of the combined data set using the WLSS 
regression yields: 

/ \ 0.21±0.03 

/!o a s (V0.7) 3/2 = (9.3^)xl0- 2 [ 2x ^: Mq 

(12) 

with a scatter of 17 per cent about the best fitting re- 
gression line. The data and resulting fit are shown in 
Figure. O As discussed in the introduction we assume 
that this relation is not evolving, in the absence of obser- 
vations to the contrary. To better elucidate trends with 
total mass, we divided the data set into the same loga- 
rithmic bins in total mass as for the stellar mass fraction, 
and computed the mean and standard deviation of the 
distribution of the gas mass fraction values in each bin. 
The large points with error bars show the trend of these 
binned data with total mass. The observed relation sug- 
gests that lower mass systems have proportionally less 
gas th an high mass system s. Further discussion is avail- 
able in lPratt et all (|2009f ). 

4.3.2. The baryon mass fraction (in galaxies and ICM) 

We now combine the results on the stellar and gas 
mass fractions derived in the previous two sections to 
investigate the behaviour of the baryonic mass fraction 
as a function of total mass. At this stage no contri- 
bution is considered from the ICL as defined in ^4.4.31 
In each logarithmic mass bin we sum the mean con- 
tribution from stellar and ICM mass components. As 
we wish to determine the behaviour of the average sys- 
tems in a given mass bin, for each component the un- 
certainty is calculated from the standard deviation of 
the mean (the standard deviation divided by yN—l, 
where N is the number of data points in the bin). The 
uncertainty on the total baryon mass content is then es- 
timated from the quadratic sum of the individual un- 
certainties for the stellar and ICM contributions. Fig- 
ure [6] (lower panel) reproduces the average behaviour of 
the sum of the two baryonic components estimated in 
the previous sections (i.e. ICM gas and stars associated 
with galaxies) as a function of total mass for galaxy sys- 
tems with 2xl0 13 <M 500 <8.1xl0 14 M©. The ensuing 
baryon mass fraction is an increasing function of the sys- 
tem mass: 

, , \ 0.09±0.03 

M500 



/500 rS+SaS = (°' 123 ± °- 003 ) X 



2 x 10 14 Mt 



(13) 

This expression is obtained after excluding the lowest 
mass point which is affected by an extremely large uncer- 
tainty since the corresponding gas fraction is estimated 
from only two groups. 

4.4. Comparison with WMAP 
4.4.1. Raw values 
As Figure [5] shows, there is a gap between the val- 
ues of /5oo r8+SaS estimated from WMAP 5 and those ob- 



tained here; this discrepancy, before any correction, is 
significant at more than 5<r for systems less massive than 
~10 14 M Q (see Table [3]), where the uncertainties are cal- 
culated as described in ^4.3.21 

4.4.2. Values corrected for gas depletion 

We now correct the value o f the baryon fracti on for 
gas depletion. As discussed in iFrenk et a l. (1999), sim- 
ulations without feedback suggest that the ICM has a 
slightly more inflated dis tribution than the dark matter 
(see also observations bv lPratt fc Arnauc3l2002f ). result- 
ing in a decrease in the gas fraction of 10% at Rsoo- 
In the absence of indications to the contrary we do not 
assume a mass dependence for the gas depletion. For 
average massive clusters ((M 500 ) = 7 x 1O 14 M ) the 
value of gas depletion-corrected ystars+gas+dcpi j s con _ 
sistent within 1.4er with the WMAP 5 estimate. However 
the gas depletion corrected value in the group regime 
((M500) = 5x 1O 13 M ) is still 4.5cr discrepant from that 
ofWMAP5 28 . 

4.4.3. Values corrected for gas depletion and ICL 

The existence of a diffuse stellar component in galaxy 
groups/clusters is now a well established observational 
result, but the way the ICL is defined and measured is 
not unique (see iZibettil 120081 for a recent review). The 
quality of our observations is insufficient to measure the 
contribution of diffuse, very low surface brightness light 
(>25.8 K-mag arcsec -2 ) within rsoo directly for indi- 
vidual systems in the sample. To quantify the amount of 
stellar mass which is associated with diffuse light that es- 
capes detection durin g the standard pho tometry extrac- 
tion with SExtr actor (|Capak et al.ll2007f ). we are guided 
by p revious observationa l results. In particu l ar we con- 
sider iZibetti et alj (120051), iKrick k. Bernstein] (|2007l ) and 
iGonzalez et al.l (|2005f) . IZibetti et all (|2005f ) used stack- 
ing analysis of 683 systems at z=0.2-0.3 ranging in total 
mass from a few times 10 13 to 5xl0 14 M Q (the average 
total mass is 7xl0 13 M Q ), selected from a 1500 deg 2 of 
SDSS-DR1, reaching the unprecedented surface bright- 
ness limit of ^32 mag arcsec -2 (R-band in the z=0.25 
observed frame). They show that on average the ICL 
contributes ~11% of th e stellar light within 500 kpc. In 
a complementary study IKrick fc Bernstein! (|2007[ ) used a 
sample of massive clusters with a range of morphology, 
rcdshift and densities to find that the ICL contributes 
with 6%-22% to the total cluster light in r-band within 
one quarter of the virial radius, finding no appreciable 
correlation with cluster mass. Given these results, we as- 
sume that the contribution of the ICL to the total mass 
of a system is equal to its observed contribution to the 
total light and ranges between 11 and 22 %. This range is 
consis tent with the theo r etical results bv lMurante et all 
(|2007h and lPurcell et alj (|2008[) . in their attempt of mod- 
elling the ICL by numerical simulations. Furthermore 
given the complete lack of observational constraints, we 
assume that the ICL mass fraction is not evolving with 
rcdshi ft for 0<z< l; this is supported b y the simulation 
of iDubinski et alj (|2003l ) as shown in iFeldmeier et al.l 

28 We note that this discrepancy represents a lower limit if a fur- 
ther 10% reducti on of the gas mas s is applied due to the dumpiness 
of the ICM as in Lin et al. (2003). However this correction is not 
applied in most of the studies of gas component in clusters. 
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TABLE 3 

Discrepancy of /b from the WMAP5 value in sigma 

UNITS 



M500/ [hra M ©1 




A/ b / [°7 6 ] a 


A/ b / k/J b 


2.1e+13 


>1.2 


>0.8 


>0.3 


5.1c+13 


5.3 


4.5 


3.3 


1.2e+14 


5.1 


4.2 


3.2 


3.0e+14 


3.7 


2.6 


2.1 


7.1o+14 


2.(5 


1.4 


1.0 



a After correction for gas depletion. 15 After correction for 
gas depletion and ICL. 
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Fig. 6. — Lower panel: average stellar to dark mass ratio (filled 
points) for the COSMOS+LM03 sample and average gas fraction 
(empty points). Uncertainties are computed from the standard de- 
viation of the mean in all cases. Upper panel: total baryonic frac- 
tion obtained summing the points in the lower panel compared with 
the universal value by WMAP5 (dark grey stripe). The dashed- 
dotted line represents the fit to the measured points. The dashed 
line represents the fit to the points where the gas fraction has been 
corrected for a 10% gas depletion. The light grey stripe is the fit 
to the relation taking in account both gas depletion and a constant 
(11-22%) ICL contribution to the stellar mass. 

( 2004). We discuss the impact of our choice on the results 
in §4.51 The final gas depletion corrected values includ- 
ing the ICL contribution of f^re+gas+depl+iCL are l owe r 
than the WMAP5 estimate across the entire explored 

f3tars+£ras+dcpl+ICL • , .m ,i 

500 is m agreement with the 

WMAP5 result within la in the massive cluster regime, 

but still discrepant at a significance level of at least 3.3c 

for groups (see Figure 0- 

4.5. Impact of systematic effects 

The basic observational result of the present study is 
that the baryon mass fraction, corrected for gas deple- 
tion and ICL contribution, is consistent with WMAP5 
estimate within la for clusters with (M)=7x 10 14 M Q 
but is significantly (3.3er) lower for groups with 
(M)=5x 1O 13 M . At the cluster scale our result on 
the ba ryon fraction is consistent with that of iLin et al.l 
(2003), indicating that different approaches do not show 



systematic differences in the determination of the gas 
fraction scaling with the cluster mass. Furth ermore we 
note t hat the scaling relation determined bv lPratt et al.l 
( 2009) is based on three different samples of groups and 
clusters: this should re duce the po t ential bias produced 
by sample selection. In lPratt et al.l (|2009f ) the best fit re- 
lation to the combined data from hydrostatic estimates 
reproduces the REXCESS sample distribution where the 
g as masses have been estimated using the M-T relation 
of lArnaud et al.l (|2005f ) . This suggests that potential sys- 
tematic effects on our estimates of the gas mass fractions 
at low redshifts are negligible. 

In the absence of direct estimates of the gas fraction 
at z>0.2, we have to rely upon the results of existing 
simulations, which predict the gas fraction within r5oo 
to increase on average by ~5% (adiabatic simulations) 
or 10-20% (simulations with cooling and star forma- 
tion) between z=0 an d z=l for groups and poor clusters 
(|Kravtsov et al.1120051 ). Applying a correction to this ef- 
fect at the median redshift of the COSMOS group sam- 
ple (by 5-10% at z=0.5), the discrepancy in the baryon 
mass fraction between groups of (M)=5x 10 13 Mq and 
WMAP5 is reduced to 3.0-2.6 er. Therefore we conclude 
that systematic underestimates of the gas fraction allevi- 
ate but do not solve the discrepancy at the group scale. 

Since inside groups the stellar mass fraction is com- 
parable to the gas mass fraction, we analyze the im- 
pact of the ICL fraction and the adopted stellar mass- 
to-light ratio (M/L) of the galaxy population. We 
have adopted a mass independent correction to the to- 
tal stellar mass fraction for ICL, equal to 11-22%. If 
a strong anti-correlation between the ICL mass fraction 
and the total mass of the system exists, and the true 
ICL mass fraction is equal to ~50% at the group scale, 
an agreement between our total baryonic mass fraction 
and the WMAP5 estimate is rea c hed. Such a figure has 
been claimed by iGonzalez et al.1 (|2007h for a sample of 
23 BCG-dominated clusters and groups. However the 
ICL-to-BCG light ratio (ICL/BCG) is strongly depen- 
dent on the decomposition of the total surface bright- 
ness profil e of the two components and the photomet- 
ric d epth (IGonzalez et al.ll2005l : IZibettil [2008) . We note 
that IGonzalez et al.l I^OOoT give ICL/BCG>5 by apply- 
ing a simultaneous decomposition of the surface bright- 
ness distribution of BCG+ICL in two De Vaucouleurs 
components: the outer one is considered as the gen- 
uine ICL and the inner one as the BCG. Conversely, 
IZibetti et all (|2005T) obtain ICL/BCG<0.5 by fitting only 
the inner profile with a De Vaucouleurs model (which 
represents the BCG) a nd consi d ering all the residual light 
as ICL. Nevertheless, IZibet ti (2008) applied a two-De 
Vaucouleurs decomposition to the IZibetti et al.l (2005) 
data obtaining ICL/BCG^2, and concluded that the 
ICL+BCG-to-total light ratio is a much more robust 
measurement, which is likely equal to 0.3 (in light) for 
systems of average mass 7xl0 13 Mf T ). Th e high value of 
ICL/BCG found bv lGonzalez et~aT1 pool may be the re- 
sult of a sample bias, as suggested by the same authors. 
On the other hand the lack of trends reported by Zibetti 
et al. (2005) could be intrinsically biased by the adoption 
of a fixed metric aperture of 500 kpc, which correspond to 
smaller fraction of R200 for more massive clusters. Given 
the steeper prole of the ICL with respect to galaxies, the 
ICL fraction of more massive clusters could be overes- 
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timated and a correction for this effect could reconcile 
these results with the negative trend found by Gonza- 
lez et al. (2007), but not with the extreme values of 
ICL+BCG-to-total light ratio. Generally it is evident 
that better determinations of the trends of the ICL with 
cluster mass and richness are needed. 

The COSMOS groups sample contains a whole range 
of systems, which exhibit a BCG-to-galaxy stellar mass 
ratio from 0.2 to 0.9. For these groups the estimated 
ICL+BCG-to-total light ratio for the average group is 
0.36, broadly consi stent with the g enerall y accep ted av- 
erage value of 0.3 ([Gonzalez et all 120071 : TZibettJ l2008T >. 
This suggests that we are not missing an important con- 
tribution of the stellar mass in our analysis, in spite of 
our definition of ICL. 

Another systematic effect may be introduced by the 
computation of the stellar mass-to-light ratio for the 
ensamble of the member galaxies and the ICL. In our 
case we use M/L values that correspond to the individ- 
ual star formation histo ries of individual member galax- 
ies (|Arnouts et al.|[2~007t ) and we do not make assumption 
on the M/L of the ICL. Hence the major source of sys- 
tematics on the stellar mass-to-light ratio of our galaxies 
is given by the adopted initial mass function (IMF). For 
instance, a change from a standard Salpeter to a Chabrier 
IMF reduces the M/L by 30% (jLonghetti k, Saraccol 
2009). This translates into a decrease by 30% of the 
stellar mass associated with galaxies which makes the 
bulk of the total stellar mass in our systems. There 
is no compelli ng reason to abandon the Salpeter IMF 
(R enzinil |2005[). but i t is a possibility explored in the 
literature. iLin et al.l ([2003) obtained the stellar mass- 
to-light ratio for the ensamble of group/cluster mem- 
ber galaxies by folding in a morphological type depen- 
dent M/ L with the temperature dependence of the spiral 
fraction; iGonzalez et all (|2007l ) assumed that the ICL 
and all member galaxies share the same stellar M/L, 
as the one that characterise an early type galaxy. The 
latter case assumes that the intergalactic stars are ho- 
mogeneous with the BCG stellar population. However, 
it has been suggested that the ICL may (also) ori- 
gin from the s tripping of non-BCG galaxies ins i de th e 
group/cluster (jPurcell et all 120081 : iPierini et all l2008h . 
which are on average bluer than the BCG, especially 
m groups ( Zabl udoff & Mulchaevj | l998l : IWeinmann et all 
120061: iPoggiantiet all l2006f).~For example, if the ICL 
mass-to-light ratio used in lGonzalez et all (|2007h is over- 
estimated by a factor 2, it translates in tghe systematic 
overestimation of the baryon mass fraction by 10%. 

This systematic effect has the same amplitude, but 
opposite sign, of the potential offset applied to the 
gas fr action-mass relation according to IGonzalez et al.l 
(2003). Therefore we conclude that a 3cr discrepancy be- 
tween the baryon mass fraction of groups and the WMAP 
value holds against major systematic effects on the stellar 
populations cither diffuse or associated with galaxies. 

An overestimate of the total M/L is not enough to ex- 
plain the values of the stellar mass fractio n for the low- 
est mass systems in IGonzalez et all (|2007h which largely 
exceed the constraint on the t otal baryon fractio n set 
by WMAP 5 (as also noted in iBalogh et~all 12008ft . A 
way out is a systematic and large underestimate of 
the total masses of these systems, as also suggested by 



IBalogh et all (|2008h 29 . 

We conclude that a robust estimate of the total mass is 
crucial for systems with the lowest mass (in our sample 
(M500) ~2x10 13 Mq). Our estimates are based on the 
LX-M200 relation established via the weak lensing anal- 
ysis in Leauthaud et al. (in preparation), and exhibit 
a typical uncertainty of 30%. The use of different total 
mass estimators could offer a test of the presence of sys- 
tematics, but unfortunately this is still hard to achieve 
for statistical large samples of groups at different red- 
shifts. 

5. DISCUSSION 

We have investigated if the discrepancy between esti- 
mates of the total baryon mass fraction obtained from 
observations of the CMB and of galaxy groups per- 
sists when a large, unbiased sample of well-characterized 
groups is considered. The COSMOS 2 deg 2 survey 
meets this requirement, yielding 91 candidate X-ray 
groups/poor clusters at redshift 0.1 < z < 1. In order to 
extend the span in total mass to two orders of magnitude 
(2xl0 13 <M 500 <1.2xl 15 M m ), w e cons ider 27 nearby 
clusters investigated by ILin et ail ((2003) . Comparable 
robust measurements of total mass and total stellar mass 
(in galaxies) exist for individual objects of both subsam- 
ples, as shown in the previous sections. In addition, the 
same scaling relation is used to estimate the gas mass 
fraction in both subsamples. This enables us to build a 
joint sample of 118 X-ray selected groups and clusters at 
z < 1 for which the importance of systematics is reduced 
(see $2]) . For this sample, the behaviour of the total stel- 
lar mass fraction as a function of the total mass can be 
investigated for a large range in total mass and, for the 
first time, in redshift (at least for groups). The results 
of our analysis and their impact on the widely accepted 
paradigm of the hierarchical growth of structure in the 
universe are discussed hereafter. 

5.1. The stellar mass fraction 

We have shown (Figure[5]) that the stellar-to-total mass 
ratio in COSMOS groups and in 27 local clusters is an- 
ticorrelated with the total mass of the system. This re- 
lation is given by /|Qo rs °< ^500' 37±a04 , which holds also 
after introducing the mass independent correction for the 
ICL (see ST}- The global trend between /|^ rs and M 500 
is consistent wi th that obsery e d in c lusters at z<0.3 both 
by LMS03 and lLagana et all ((2008) using much smaller 
samples. We extend their results to the low mass regime 
by one decade and to higher redshift. 

The difference in the number of stars formed per unit 
of halo mass between groups and clusters has been in- 
terpreted in terms of a varying efficiency of th e star for- 
matio n with the total mass of the system (e.g. ILin et all 
2003). A variation in the star-formation efficiency for 
systems with v irial temperatures > 10 7 K is a result of 
simulations bv ISpringel fc HernquistJ (|2003h : it is inter- 
preted in terms of cooling flows being less efficient in 
shutting off star formation in groups. An alternative 
possibility is that clusters are formed not only by merg- 
ing of groups and smaller clusters but also that they ac- 

29 These objects certainly impact the strongly inverse total 

mass dependen ce of the total stellar mass fraction found by 

IGonzalez et~aTl [12001 ). 
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crete a large fraction of their galaxies (with a low stel- 
lar mass fraction, of the order of 0.01) from t he field 
dWhite fc Frenld[l99ll : [Marinoni fc Hudson! 12002ft . How- 
ever after a mass independent correction for the ICL 
contribution (introduced in §4.4j) . the relation /foo rs oc 
MgQQ 37 * ' 04 is in agreement with the constraint on the 
slope set by the hierarchical model of structure forma- 
tion under the assumption that a t least half of the st ars 
in groups were formed by z = 1 (|Balogh et al.l [20 08) 30 . 
This is supported by the apparent absence of evolution 
for this relation in our sample within the redshift range 
0.1-1. This shows how observational studies such as the 
present one can improve the constraints on models and 
foster our understanding of the underlying physical pro- 
cesses. 

5.2. The total baryon mass fraction 

Combining the computed stellar mass fraction with the 
estimated gas mass fraction deri ved from the mean lo- 
cal relation in iPratt et al.l (|2009f l. we find that the gas 
plus stellar (galaxies) baryon mass fraction increases by 
^25% (from ~0. 1 1 to ^0.14) when the total mass in- 
creases by a factor of one hundred. After a constant 
10% correction for gas depletion and a further correc- 
tion for a constant 11-22% ICL contribution, the value 
°f /g*^ IS+gas+dcpl+ICL for an average cluster is consistent 
within lcr with the cosmic value measured by WMAP, 
while the jstars+gas,+dcpi+iCL found for an average group 
differs from it at more than 3a. Given the heterogeneity 
of the sample (see e.g. Figure [5]), for some objects the 
gap between ^tare+gas+depl+iCL and ^ WMAp5 value 

could be negligible or, conversely, statistically more sig- 
nificant for objects in the same bin of total mass, but 
at the two extremes of the distribution in fstms- Un- 
fortunately we do not have a measure of the gas mass 
fraction for individual objects, therefore we focus on the 
behaviour of the average object. We did likewise for the 
ICL by assuming a fixed fractional contribution of 11- 
22% across the entire mass range. Possible systematic 
effects introduced by our definition and estimate of the 
ICL contribution are discussed in ij4.5l Here we stress 
that they do not lead to an anomalously low BCG+ICL 
contribution to the total mass of the system. Thus the 
discrepancy at the groups regime in not erased by uncer- 
tainties on the stellar mass fraction. In the absence of 
evidence for a systematic and relevant underestimation 
of the gas mass fraction in our systems (see ^4.5|) , we in- 
terpret the discrepancy as a lack of gas, by 33%, at the 
group regime. This may be produced by feedback (stellar 
and/or AGN), as suggested by high-resolution cosmolog- 
ical simulations including cooling, star formation, super- 
nova feedback, and A GN radio-mode feedback in galaxy 
clust e rs and groups dPuchwein et al.ll2008l iBower et all 
120081 iShort & Tholnasl feoOS) . Since supernova feedback 
appears to be insuffici ent to explain the Lx~T relation 
(|Puchwein et al.|[2008h . feedback by AGN seems neces- 
sary. According to this interpretation, gas can be re- 
moved from within R500 mainly as a consequence of the 
mechanical heating produced by a central AGN. The ac- 
tion of the AGN is larger in groups than in clusters sim- 

30 We note that a steeper relation is o btained when the stron gly 
inverse mass dependen t ICL fraction of I Gonzalez et al. (2007) is 
used (sec Balogh ct al. (2008) for the discussion). 



ply because the potential well is shallower in the former 
systems. In a forthcoming work we will quantify the 
feedback by AGN radio-mode for the COSMOS groups. 
Another proposed mechanism capable of accounting for 
the " missing" gas is "filamentary heating" (|Voit fe; Brvanl 
2001). Low entropy gas is consumed in star formation 
before the group formation, which eventually raises the 
entropy of the gas which becomes the ICM. The resulting 
higher entropy level inhibits the gas from falling towards 
the center of the potential well, which can e xplain the 
lack of gas in the central region of groups (Su n et aT] 
120091) . 

6. CONCLUSIONS 

The baryon mass fraction is a parameter which can be 
constrained by the primordial light elements abundance 
set by the nucleosynthesis at early epochs. It can be inde- 
pendently measured from observations of the CMB (e.g. 
WMAP) or of galaxy groups/clusters. Different studies 
of the baryon mass fraction in nearby galaxy systems 
have reported values lower than the one from WMAP, 
the discrepancy being larger for groups than clusters. We 
investigate if this discrepancy persists when a sample of 
local clusters is supplemented by a large, unbiased sam- 
ple of groups at 0.1 < z < 1.0. Hereafter we list our 
conclusions. 

1. The stellar mass fraction associated with galax- 
ies is anticorrelated with the mass of the system: 
/f£ff s cx M^ ° 37±0M . This is consistent with pre- 
vious results on local clusters. The validity of this 
result is now extended by one decade in total mass 
and to redshift 1. 

2. The previous relation holds after correcting the 
stellar mass fraction for a mass independent 11- 
22% contribution from the ICL as suggested by 
both observations and simulations. The slope of 
the /|Jjo rs -Af5oo relation is consistent with the con- 
straint set by t he hierarchical par adigm of struc- 
ture formation (|Balogh et al.ll2008h . No significant 
evolution in the relation between f|go rs and M 50 o is 
observed. This supports the scenario in which mas- 
sive clusters form mostly by merging of less massive 
groups and clusters, and observed groups in the 
redshift range 0-1 have formed the bulk of their 
stellar mass by z^l.0. 

3. Combining measured values of the stellar mass 
fraction with values of the gas mass fraction esti- 
mated from an average relation obtained for a local 
sample, f^j^+sas increases by 25% from groups to 
clusters. After the introduction of appropriate cor- 
rections for gas depletion and ICL contribution, the 
total baryonic mass fraction at the groups regime 
still differs from the WMAP 5 value at 3.3er. We in- 
terpret the origin of this discrepancy as a lack of gas 
(by 33%), which can be produced either by feed- 
back (supernovae and/or radio-mode AGN heat- 
ing) or by "filamentary heating". 

Our results provide useful constraints on simulations of 
the aforementioned processes. In particular the availabil- 
ity of a large unbiased sample of groups offers direct and 
stringent constraints on models rather than relying on 
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extrapolation of the behaviour of the stellar fraction as 
a function of mass in the entire family of systems with 
10 13 <M 50 o < 10 14 M Q . Future observations will increase 
both the statistics and the redshift sampling rate, so that 
a test and extension of our conclusions will be possible. 

The authors thank the anonymous referee for her/his 
valuable comments, which led to a significant im- 
provement of the paper. We acknowledge the 



contributions of the entire COSMOS collaboration; 
more informat ions on the COSMO S survey are avail- 
able at |http : //www . astr . caltech . edu/ $\sim$cosmos. 
This research was supported By the DFG Cluster 
of Excellence Origin and Structure of the Universe 
(http://www.universe-cluster.de). D.P. acknowl- 
edges support by the German Deutsches Zentrum fur 
Luft- und Raumfahrt, DLR project number 50 OR 0405. 



REFERENCES 



Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2001, MNRAS, 328, 
L37 

Allen, S. W., Schmidt, R. W., Ebeling, H., Fabian, A. C, & van 

Speybroeck, L. 2004, MNRAS, 353, 457 
Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 

893 

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, 
L37 

Arnouts, S., et al. 2002, MNRAS, 329, 355 
Arnouts, S., et al. 2007, A&A, 476, 137 
Bahcall, N. A. 1977, ARA&A, 15, 505 

Balogh, M. L., McCarthy, I. G., Bower, R. G., & Eke, V. R. 2008, 

MNRAS, 385, 1003 
Bartelmann, M. 1996, A&A, 313, 697 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 

Bialck, J. J., Evrard, A. E., & Mohr, J. J. 2001, ApJ, 555, 597 

Biviano, A., & Salucci, P. 2006, A&A, 452, 75 

Bower, R. G., Benson, et al. 2006, MNRAS, 370, 645 

Bower, R. G., McCarthy, I. G., & Benson, A. J. 2008, MNRAS, 

390, 1399 
Brown, J. P. 1997, Ph.D. Thesis 

Bundy, K., Ellis, R. S., & Conselice, C. J. 2005, ApJ, 625, 621 
Capak, P., et al. 2007, ApJS, 172, 99 

David, L. P., Arnaud, K. A., Forman, W., & Jones, C. 1990, ApJ, 
356, 32 

Diehl, S., & Statler, T. S. 2007, ApJ, 668, 150 

Dressier, A. 1980, ApJS, 42, 565 

Dressier, A., & Schectman, S. A. 1988, AJ, 95, 284 

Dubinski, J., Koranyi, D., & Geller, M. 2003, Astrophysical 

Supercomputing using Particle Simulations, 208, 237 
Dunkley, J., et al. 2009, ApJS, 180, 306 

Elvis, M., & C-COSMOS Team 2006, Bulletin of the American 

Astronomical Society, 38, 1006 
Ettori, S. 2003, MNRAS, 344, L13 
Evrard, A. E. 1997, MNRAS, 292, 289 

Feldmeier, J. J., Mihos, J. C, Morrison, H. L., Harding, P., Kaib, 

N., & Dubinski, J. 2004, ApJ, 609, 617 
Finoguenov, A., Burkert, A., & Bohringer, H. 2003, ApJ, 594, 136 
Finoguenov, A., et al. 2007, ApJS, 172, 182 
Frenk, C. S., et al. 1999, ApJ, 525, 554 

Gilfanov, M., Grimm, H.-J., & Sunyaev, R. 2004, MNRAS, 351, 
1365 

Gonzalez, A. H., Zabludoff, A. I., & Zaritsky, D. 2005, ApJ, 618, 
195 

Gonzalez, A. H., Zaritsky, D., & Zabludoff, A. I. 2007, ApJ, 666, 
147 

Grimm, H.-J., Gilfanov, M., & Sunyaev, R. 2003, MNRAS, 339, 
793 

Guzzo, L., et al. 2007, ApJS, 172, 254 

Hasinger, G., et al. 2007, ApJS, 172, 29 

He, P., Feng, L.-L., Fang, L.-Z. 2006, ApJ, 623, 601 

Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 

Ilbert, O., et al. 2009, ApJ, 690, 1236 

Koekemoer, A. M., et al. 2007, ApJS, 172, 196 

Komatsu, E., et al. 2009, ApJS, 180, 330 

Kravtsov, A. V., Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588 

Krick, J. E., & Bernstein, R. A. 2007, AJ, 134, 466 

Kiipcii Yoldas,, A., et al. 2007, A&A, 463, 893 

Lagana, T. F., Lima Neto, G. B., Andrade-Santos, F., & Cypriano, 

E. S. 2008, A&A, 485, 633 
Leauthaud, A., et al. 2007, ApJS, 172, 219 
Lilly, S. J., et al., 2007, ApJS, 172, 70 

Lin, Y.-T., Mohr, J. J., & Stanford, S. A. 2003, ApJ, 591, 749 



Lin, Y.-T., & Mohr, J. J. 2004, ApJ, 617, 879 
Longhetti, M., & Saracco, P. 2009, MNRAS, 174 
Marinoni, C, & Hudson, M. J. 2002, ApJ, 569, 101 
Massey, R., et al. 2007, Nature, 445, 286 

McCarthy, I. G., Bower, R. G., & Balogh, M. L. 2007, MNRAS, 
377, 1457 

Murante, G., Giovalli, M., Gerhard, O., Arnaboldi, M., Borgani, 

S., & Dolag, K. 2007, MNRAS, 377, 2 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1995, MNRAS, 

275, 720 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 

493 (NFW) 
Pannella, M., et al. 2006, ApJ, 639, LI 

Pierini, D., Zibetti, S., Braglia, F., Bohringer, H., Finoguenov, A., 

Lynam, P. D., & Zhang, Y.-Y. 2008, A&A, 483, 727 
Poggianti, B. M., et al. 2006, ApJ, 642, 188 
Pratt, G. W., & Arnaud, M. 2002, A&A, 394, 375 
Pratt, G. W., Arnaud, M., & Pointecouteau, E. 2006, A&A, 446, 
429 

Pratt, G. W., Croston, J. H., Arnaud, M., & Bohringer, H. 2009, 
A&A, 498, 361 

Puchwein, E., Sijacki, D., & Springel, V. 2008, ApJ, 687, L53 
Purcell, C. W., Bullock, J. S., & Zentner, A. R. 2008, MNRAS, 
391, 550 

Reiprich, T. H., & Bohringer, H. 2002, ApJ, 567, 716 
Renzini, A. 2005, The Initial Mass Function 50 Years Later, 327, 
221 

Roussel, H., sadat, R., & Blanchard, A. 2000, A&A, 361, 429 
Rykoff, E. S., et al. 2008, MNRAS, 387, L28 
Sanders, D. B., et al. 2007, ApJS, 172, 86 

Sales, L. V., Navarro, J. F., Lambas, D. C, White, S. D. M., & 

Croton, D. J. 2007, MNRAS, 382, 1901 
Salpeter, E. E. 1955, ApJ, 121, 161 
Salvato, M., et al. 2009, ApJ, 690, 1250 
Sandage, A. 1986, A&A, 161, 89 
Schechter, P. 1976, ApJ, 203, 297 
Schombert, J. M. 1986, ApJS, 60, 603 
Scoville, N, et al. 2007, ApJS, 172, 1 

Scoville, N, et al. 2007, ApJS, 172, 150 

Short, C. J., & Thomas, P. A. 2008, arXiv:0811. 3166 

Springel, V., & Hernquist, L. 2003, Astrophysical Supercomputing 

using Particle Simulations, 208, 273 
Sun, M., Voit, G. M., Donahue, M., Jones, C, Forman, W., & 

Vikhlinin, A. 2009, ApJ, 693, 1142 
Taniguchi, Y., et al. 2007, ApJS, 172, 9 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, 

M., Murray, S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 
Voit, G. M., & Bryan, G. L. 2001, Nature, 414, 425 
Weinmann, S. M., van den Bosch, F. C, Yang, X., & Mo, H. J. 

2006, MNRAS, 366, 2 
White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 
White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 

1993, Nature, 366, 429 
Wilkins, S. M., Trentham, N., & Hopkins, A. M. 2008, MNRAS, 

385, 687 

Zabludoff, A. I., & Mulchaey, J. S. 1998, ApJ, 496, 39 
Zamojski, M. A., et al. 2007, ApJS, 172, 468 

Zhang, Y.-Y., Finoguenov, A., Bohringer, H., Kneib, J. -P., Smith, 
G. P., Kneissl, R., Okabe, N, & Dahle, H. 2008, A&A, 482, 451 

Z hao, D. H., Jin g, Y. P., Mo, H. J., & Boerner, G. 2008, 
arXiv:0811.0828 

Zibetti, S. 2008, in Dark Galaxies and Lost Baryons, IAU 
Symposium, Vol. 244, p. 176-185 



14 



Giodini ct al. 



Zibetti, S., White, S. D. M., Schneider, D. P., & Brinkmann, J. 
2005, MNRAS, 358, 949 



